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We investigate numerically kinematic dynamos driven by flow of electrically conducting fluid in the shell 
between two concentric differentially rotating spheres, a configuration normally referred to as spherical 
Couette flow. We compare between axisymmetric (2D) and fully three dimensional flows, between low and 
high global rotation rates, between prograde and retrograde differential rotations, between weak and strong 
nonlinear inertial forces, between insulating and conducting boundaries, and between two aspect ratios. 
The main results are as follows. Azimuthally drifting Rossby waves arising from the destabilisation of the 
Stewartson shear layer are crucial to dynamo action. Differential rotation and helical Rossby waves combine 
to contribute to the spherical Couette dynamo. At a slow global rotation rate, the direction of differential 
rotation plays an important role in the dynamo because of different patterns of Rossby waves in prograde 
and retrograde flows. At a rapid global rotation rate, stronger flow supercriticality (namely the difference 
between the differential rotation rate of the flow and its critical value for the onset of nonaxisymmetric 
instability) facilitates the onset of dynamo action. A conducting magnetic boundary condition and a larger 
aspect ratio both favour dynamo action. 
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Introduction 



Dynamo action is believe d to generate magnetic fields in the universe (e.g. 
Riidiger and Hollerbach 20051 ). In dynamo action, the motion of conducting fluid produces 
a magnetic field through the effect of electromagnetic induction while the generated magnetic 
field gives a back-reaction on the fluid motion, and so the dynamo is a coupled nonlinear sys- 
tem. A simplified model is the kinematic dynamo in w hich the back-reaction of the field on the 
flow is neglected. Beginning with the seminal work of Bullard and Gellmanl (| 19541 ). 

kinematic 

dynamo theory is a classic test-bed for discovering the efficien cy of flows in generating mag- 



netic field. The subject has already been extensively stu died (jRobertsI 119721 . iGubbins et al 



2000al lbh and is most recently revi ewed bvlGubbind (l2008h. In the kinematic dynamo the fluid 
flow can either be prescribed (e.g. Galloway and Proctor! 19921 ) or be a solution to the Navier - 
Stokes (N-S) equation of fluid motion subject to a given forcing (e.g. Livermore et al. 20071 ). 
In the latter case of prescribed forcing the nonlinear interaction between flow and field is re- 
moved by dropping the Lorentz force in the N-S equation. Because in the kinematic dynamo 
the magnetic field has no back-reaction on the fluid flow, the magnetic induction equation is 
linear and thus the magnetic field does not saturate but grows. 

In the geodynamo driven by convective fluid motion in the Earth's core, it is plausible 
that a large scale differential rot ation in the zonal flow emerges because of angular momen- 
tum transfer through convection (|Bussdll978l . iManneville and Qlsonl ll996. lAurnou and Olsonl 
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1200 ih . i.e. the so-called thermal wind, or through turbulent Reynolds stress (jRiidigerlll989l ) 



such that the solid inner core rotates at a different rate relative to the mantle due to the viscous 
coupling exerted on the inner core. This large scale differential rotati on has been supported 
by th e recent geodynamo simulations at a very low Ekman number ( Sakuraba and Roberts! 
20091 ) and the super-rotation of the s olid inner core by some seism ic data ( Zhang et al. 20051 ) 
(though it still remains controversial (jSouriau and Poupinetll2003l )). This differential rotation 
is important in dynamo action, e.g. through the u effect that shears poloidal magnetic field 
lines to create toroidal magnetic field lines. Many simulations of the geodynamo driven by 
thermal and compositional convection have been carried out in the geomet ry of a spheri- 
cal shell ( Roberts and Glatzmaierl 2000l . Wicht and Tilgner 2010l . Jones! 2011 ). However, the 
convectively driven dynamo is very complicated, and therefore a simplified dynamo model is 
proposed, i.e. the spherical Couette dynamo, in which the dynamo is driven by an imposed 
mechanical force arising from the differential rotation of the inner and outer spheres. This 
simple dynamo model captures one of the important ingredients in dynamo action, namely 
differential rotation. The need to understand dynamo action at realistic values of the magnetic 
Prandtl number has spawned a number of liquid metal dynamo exper iments using differentia l 
rotation in a sp herical geometry, e xamples of which are in Grenoble ( Cardin and Bri to 2007) 
and Maryland ( Pefflev et al. 2000). These experiments motivate the present study. 

Flow in a spherical shell driven by differentially rotating inner and outer spheres is called 
spherical Couette flow. Sph erical Couette flow at an infinitesimally differential rotation rate 
was analytically studied by Stewartson ( 19661 ). He showed that the fluid outside the tangent 
cylinder, a cylinder parallel to the rotational axis and touching the inner sphere, rotates at the 
rate of the outer sphere due to it experiencing a homogeneous boundary condition, whereas 
the fluid inside the tangent cylinder rotates at an intermediate rate between the rotation rates 
of inner and outer spheres due to it experiencing an inhomogeneous boundary condition. Con- 
sequently a shear layer in the angular velocity emerges at the tangent cylinder, which is called 
the Stewartson layer. Subsequently, the situat ion of finite differ ential rotation rates was nu- 
merically studied by Hollerbach and coworkers. Hollerbach! ( 20031 ) calculated the axisymmetric 
basic state and the nonaxisymmetric linear stability of spherical Couette flow, and pointed out 
that the instability of Stewartson flow is asymmetric in th e direction of differentia l rotation, 
namely retrograde flow is more stable than prograde flow. Hollerbach et al. ( 20041 ) (see also 
Schaeffer and CardiiJlioOsh calculated the nonlinear saturation of instability, and showed that 
the destabilisation of the Stewartson layer triggers azimuthall y drifting Rossb y waves due to 
the curvature of the spherical geometry, the so- called (3 effect (|Pedloskvill987h . 

In addition to nonmagnetic spherical Couette flow, the electrically conducting spherical Cou- 
ette flow with an imposed field is called the magnetic spherical Couette flow . The axisymmetri c 
basic state of m agnetic spher i cal Couette flow w as studied analytica ll y bvlStarchenkol (1998), 
numer ic ally by Dormv et al. ( 1998 . 20021 ) and Hollerbach ( 2000a ). IHollerbach and Skinner! 
(|200ll ). IHollerbach et al\ (I2007T), and experime ntally in Maryland (jSisan et al\ |2004 ) and 
Grenoble ( Schmitt et al. 20081 . Nataf et~a/1l2008l ) . When the imposed field is uniformly aligned 
with the rotation axis it enhances the shear in the Stewartson layer, however, when the imposed 
field is not parallel to the rotation axis (say, dipolar or quadrupolar) the increasing imposed 
field alters the Stewartson layer until a new curved shear layer emerges, called the Shercliff 
layer. The nonaxisymmetric instability in magnetic s pherical Couette flow was then stud- 
ied w ith an imposed field parallel to the rotation axis ( Wei and Hollerbach 20081 . Hollerbach 



20091 ) . It is found that the imposed field can suppress the asymmetry of instabilities between 
prograde and retrograde flows in the nonmagnetic Stewartson layer. 

We now come back to the Couette dynamo problem. Dynamo action induced by Taylor- 
Couette flow in a cylindrical geometry has been studied by Willis and Barenghi (2002) in both 
linear and nonlinear regimes, in which shear and Taylor vortices are combined to contribute 
to dynamo action. However, the two geometries are quite different; the cylindrical problem 
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has neither Stewartson layers nor a f3 effect to trigger Rossby waves, whereas the spherical 
p roblem has no Taylor vortice s . 

Schaeffer and Cardin ( 20061 ) studied the spherical kinematic dynamo in a particular form of 



rapidly rotating flow. There is no inner sphere in their spherical geometry but the outer sphere 
is split into two parts along the tangent cylinder, the region outside the tangent cylinder ro- 
tating at the global rate and the top and bottom boundaries inside rotating at a different rate. 
In this configuration the Stewartson shear layer also emerges because of the angular velocity 
jump and the destabilisation of the shear layer triggers Rossby waves. A quasi-geostrophic 
model for the flow was used, in which the cylindrical radial and azimuthal velocities are in- 
variant along the rotation axis {z direction) but the axial velocity is linearly dependent on 
the z coordinate, such that a regime of extremely low viscosity is numerically achievable (very 
low Ekman and moderately small magnetic Prandtl numbers). In purely hydrodynamic cal- 
culations they found that for prograde flow Rossby waves tend to spread out of the shear 
layer whereas those for retrograde flow tend to be confined in the vicinity of the shear layer. 
In the kinematic dynamo calculations they found that retrograde flow slightly reduces the 
dynamo threshold. The spatial structure o f Rossby wav e s is helical an d it is well known that 
a helical wave can induce the alpha effect (Parker 19551 . Moffattl 19781 ). the effect that twists 



toroidal magnetic field lines to create poloidal magnetic field lines. Therefore this dynamo was 
interpreted with the a-oo mechanism in which the u effect is induced by the shear and the a 
effect by Rossby waves. Moreover, in their calculations both dipolar and quadrupolar fields 

w ere generated. 

Guervillv and Cardin ( 20ld ) then studied the nonlinear spherical Couette dynamo in which 



the nonlinear interaction of fluid flow and magnetic field is accounted for by retaining the 
Lorentz force in the N-S equation. They calculated both the case of a stationary outer sphere 
without global rotation and the case of a rotating outer sphere with global rotation. They 
found that the resulting axisymmetric flow cannot generate a dynamo until nonaxisymmetric 
hydrodynamical instabilities are excited, and that the presence of global rotation facilitates 
the onset of dynamo action. Without global rotation the critical magnetic Prandtl number 
Pm is of the order of unity, which may be translated into a magnetic Reynolds number Rm 
to be of the order of a few thousand. However, with global rotation the critical Pm is reduced 
by half, and retrograde flow (negative Rossby number) is better for dynamo action. More 
interestingly, a dynamo window in which Rm ~ 300 appears in the regime of high Ekman 
number E = 10~ 3 and negative Rossby number —2 < Ro < —1.5. This dynamo window was 
physically interpreted as the enhancement of dynamo action by the shear layer. 

Since the nonlinear spherical Couette dynamo is complicated to understand, the kinematic 
spherical Couette dynamo is numerically investigated in this paper, in which the interaction 
of fluid flow and magnetic field is decoupled by dropping the Lorentz force. We do not study 
the case of a stationary outer sphere without global rotation but instead focus on the case 
of a rotating outer sphere with global rotation which is relevant to the Earth. Although this 
kinematic dynamo is simple, it consists of the two most important ingredients in dynamo 
action, namely the u effect induced by differential rotation and the a effect induced by helical 
Rossby waves. In what follows we describe our numerical calculations of kinematic dynamos 
in spherical Couette flow. In section 2 we show the governing equations and introduce the 
numerical method, in section 3 we discuss the results in detail, and in section 4 we draw some 
conclusions. 
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2. Governing equations and numerical method 

Following Hollerbach ( 20031 ). the dimensionless Navier-Stokes equation of fluid motion is 



^- + \Ro\U- VU = -Vp + EV 2 U + 2U x e z , (1) 
at 

and the dimensionless magnetic induction equation is 

dB 

{Ro^ 1 -^- = V x (U x B) + Rm~ 1 V 2 B. (2) 

In (1) and (2), length is normalised with the spherical gap L = r Q — rj, time with the inverse 
of global rotation rate Q~ l , velocity with |Z\J?|L where AQ is the differential rotation rate, 
pressure with pQ\AQ\L 2 where p is the fluid density, and magnetic field with w ppfl\Afi\L 
where p is the vacuum magnetic permeability. It should be noted that the absolute value of 
differential rotation \AQ\ but not AQ itself is used in the normalisation. The reason for this 
is to ensure that the dimensional and dimensionless velocities are always parallel rather than 
anti-parallel. 

In the governing equations there are three dimensionless parameters, E, Ro and Rm. The 
Ekman number 

E Q(r Q - Ti) 2 ^ 
measures the global rotation, where v is the fluid viscosity. The Rossby number 

no -™ ( 3b) 

measures the differential rotation. The magnetic Reynolds number 

J] 

measures the electromagnetic induction effect, where rj is the magnetic diffusivity. Moreover, 
the Reynolds number Re = \ AQ\L 2 /u and the magnetic Prandtl number Pm = v/rj can be 
further deduced using these three numbers, 

\Ro\ Rm RmE 

Re = -ET> Pm = -R7 = W (3d ' 6) 

In addition to glocal Reynolds number Re we can further define Reynolds number of u effect 
Re u and Reynolds number of a effect Re a . Re u should be defined with the axisymmetric 
azimuthal flow, and because in the spherical Couette flow the axisymmetric azimuthal energy 
is nearly the total axisymmetric energy we define 



Re„ = yj ^Re, (4a) 

where EK is the total kinetic energy, EKq is the axisymmetric kinetic energy. Re a should be 
defined with the nonaxisymmtric kinetic energy 



Re a = ]j^Re, (4b) 

where EK\ is the nonaxisymmetric kinetic energy. 

To achieve a Couette dynamo in laboratory experiments, Rm should be above its critical 
value (~ 10 2 ) and Pm in many conducting liquids is very small (~ 10~ 6 ), such that either 
\Ro\ should be large enough or E should be small enough (equation ([3^)). Therefore both the 
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global and differential rotation rates as well as the length scale of fluid motion should be large 
enough, which requires strong mechanical driving. 

For comparison be tween our kinematic spheric al Couette dynamo and the nonlinear spheri- 
cal Couette dynamo (IGuervillv and Cardinll201(]h . we list below the translations of parameters 



between our and their definitions, 



E = - ~— ~t7t , Ro = — — , Rm = ^—Rm gc , (5a,b,c) 

(i-xr x x 

where the subscript "gc" denotes the definition in IGuervillv and Cardin (j2010h and x is the 
aspect ratio 0.35 in IGuervillv and Cardinl l|20lrt ). 



The velocity boundary conditions are standard no-slip conditions in Couette flow, 

Ro 

U = at r = r Q , U = r sin 9 e/> at r = rj, (6a,b) 

\Ro\ 

where 4> is the azimuthal unit vector. In most of our numerical calculations r Q and are set 
to be 1.5 and 0.5 in order to simulate the aspect ratio 1/3 of the Earth's core, except in one 
calculation in which r Q = 2 and rj = 1 are chosen to have a different aspect ratio 1/2 for 
comparison with 1/3. The magnetic boundary conditions in most of our numerical calcula- 
tions are those appropriate to matching the field to an insulating exterior (both outside the 
outer sphere and inside the inner sphere), except in one calculation in which perfectly elec- 
trically conducting boundaries are employed for comparison with the insulating boundaries. 
The magnetic boundary conditions are then 

V x B = for i, r ■ B = and fx(VxB) = for c, (7a,b) 

where r is the radial unit vector, and "i" denotes insulating boundary and "c" perfectly 
conducting boundary. 

It should be also noted that although (1) and (2) are decoupled, the solution to (1) consists of 
time-dependent azimuthally drifting waves when \Ro\ is above its critical value, and therefore 
in the calculation of 3D flow-driven dynamos we need to solve (1) a nd (2) simu l taneou sly. 

The calculations are carried out with the pseudo-spectral code of iHollerbach ( 2000bl ) . The 



toroidal-poloidal decomposition method is employed such that the divergence free condition 
of the incompressible fluid flow and of the magnetic field is automatically satisfied, 

U = V x (er) + V x V x (fr), B = V X (gf) + V X V X (hr ). (8a,b) 

The geometry is a spherical shell such that spherical harmonics are used. Take the toroidal 
scalar of velocity e for example, 



e (r, 9, cf>,t)=J2 eim (r, t) p| H (cos 9)e im t (9) 

Lm 



where pj m ' is the associated Legendre polynomials of the degree I and the order m. The 
Chebyshev polynomials are used in the radial direction, 

&im (r, t) = y2 e ^im (t) Tfc(ff), x = (2r-r - ri)/(r -n), (10a,b) 

k 

where Tfc(x) is the Chebyshev polynomials. A 2nd order Runge-Kutta method is employed 
for time integration. In our calculations, Chebyshev polynomials with truncations as high as 
170 are used, and a degree and an order of spherical harmonics as high as, respectively, 170 
and 20, are used. To search for the critical Rm for dynamo action, we vary Rm in steps of 50, 
such that the error in the determination of the critical Rm is within 50, which is sufficiently 
small compared to the typical value of critical Rm = O(10 3 ) itself. For example, in the case 
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E = 1CT 3 - 5 and Ro = -1, at Rm = 2600 the magnetic energy decays whereas at Rm = 2650 
it grows, so we know that the critical Rm is between 2600 and 2650. 

In the magnetic induction equation ([2]), the advection time scale for spherical Couette flow 
is of the order of unity and short, and therefore the time for the growth/decay rate of magnetic 
field to settle in is determined by magnetic diffusion time scale or Rm / Ro. Rm and Ro of the 
dynamo solutions are of the order of, respectively, 10 3 and unity, such that it takes around 
one month to obtain one solution with the serial code which was used in our calculations. 



Results and discussions 



In figure [T] we r eproduce the results for linear stability calculations given in figure 4 of 
HollerbachT(|2003h . showing the neutral stability curves for spherical Couette flow. The Ek- 
man and Rossby numbers defined in figure Q] are identical to our definitions in equations ((3^) 
and (J3h>)- In our calculations, two values of Ekman number, 10 -3 ' 5 and 10 -4 , are chosen for 
comparison between low and high global rotation rates. At E = 10 -3 ' 5 the critical Ro for 
the onset of Rossby waves is Ro c = +0.5 for prograde flow and Ro c = —0.8 for retrograde 
flow, while at E = 10~ 4 the critical Ro is Ro c = +0.2 for prograde flow and Ro c = —0.4 for 
retrograde flow. We firstly try to use axisymmetric flows at E = 10 -35 and E = 10 -4 to pro- 
duce a dynamo, namely Ro numbers used in these calculations are below the neutral stability 
curves in figure [1] (|i?o| < |i?o c |) such that no nonaxisymmetric instability occurs and flow is 
2D. But unfortunately, regardless of whether we use a prograde or retrograde flow, all these 
axisymmetric flows fail to produce a dynamo (in the calculations we increase Rm to 10 4 but 
no dynamo occurs). So axisymmetric flow at this Ekman number regim e is not favourable for 
dyna mo action. This is the same as in the nonlinear Couette dynamo (jGuervillv and Cardin 
2010h . 





Figure 1. The neutral stability curve of spherical Couette flow (figure 4 in (Hollerbach 2003), courtesy of Journal of 
Fluid Mechanics (Cambridge University Press)). E = u/(f2(r — r^) 2 ), Ro = AQ/Q. Left panel is for prograde flow 
(Ro > 0) and right panel retrograde flow (Ro < 0). Integers beside the curves indicate the most unstable azimuthal 
modes. Letters a, b, c and d represent the first five solutions in Table [T] 

We then increase Ro slightly above its critical value to invoke the nonaxisymmetric compo- 
nents of flow which a re azimuthally drifting; Rossby waves arising from the destabilisation of 
the Stewartson layer ( Hollerbach et adl2004l ). and we achieve some dynamos. In what follows 
we show six illustrative dynamo solutions. Table 1 gives the parameters of the six dynamo solu- 
tions shown in the figures. The critical Rm of these dynamo solutions is Rm c = O(10 3 ) and the 
magnitude of the Rossby number is \Ro\ = 0(1) (see Tabled]). We ca lculate until \Ro\ < 1; 
this is a little more modest than the nonlinear dynamo calculations (jGuervilly and Cardin 
2010h in which Ro translated to our definition reaches +8.6 for prograde flow and —11.4 for 
retrograde flow. According to equation (J3k), Pm in our calculations is between O(0.1) and 
0(1), which is usually used in most simulations of convectively driven dynamo. In the first 
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four solutions we employ insulating boundary conditions and an aspect ratio rj/r G = 1/3 (let- 
ters a, b, c and d correspond to the circles in figured]). In the fifth solution we use perfectly 
electrically conducting boundary conditions and in the sixth rj/r = 1/2. 



Figure 
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Table i. The parameters of six dynamo solutions. Column 1 gives the number of the figure 
illustrating the solutions, and letters a, b, c, and d correspond to the circles in figure [TJ Rm c 
in column 4 is an upper bound of the critical Rm. Column 6 shows the magnetic boundary 
condition, "i" denoting insulating boundary and "c" perfectly conducting boundary. Columns 
7 and 8 show the percentage (%) of, respectively, toroidal kinetic energy and toroidal magnetic 
energy. Columns 9 and 10 show the Reynolds numbers of, respectively, u effect and a effect. 
Column 11 shows the parity of magnetic field, 'D' denoting dipolar equatorial symmetry and 
'Q' quadrupolar equatorial symmetry. Columns 12, 13 and 14 show the spectral truncations 
of, respectively, Chebyshev polynomials, degree and order of spherical harmonics, number on 
the left of slash for flow U and number on the right for field B. 

The fact that a 2D flow at the Ekman number regime in our calculations cannot produce 
a dynamo but a 3D flow can indicates that azimuthally drifting Rossby waves are crucial in 
this configuration for dynamo action bec ause Rossby waves with a helical spatial structure 
can induce an a effect of dynamo action ( Avalos-Zuniga et al. 2009), namely the motion of 



helical waves twists toroidal magnetic field lines to create poloidal magnetic field lines. Here 
we need to emphasise that our usage of the term a effect is a qualitative description of the 
helical nature of the flow, but is neither connected to a kinematically described a effect nor 
to an electromotive force that results from an average over small scale correlat i ons i n the 
velocity and magnetic fields, namely the mean field theory ( Krause and Radler 1980l ). On 



the other hand, in all the dynamo solutions the toroidal kinetic energy of the zonal flow 
which is in a differential rotation occupies around 90% of the total kinetic energy, and the 
toroidal magnetic energy occupies around 80% or even higher of the total magnetic energy 
(see Table [1]) . Therefore it is reasonable to conclude that these spherical Couette dynamos 
are a-uj dynamos in which the uj effect is induced by differential rotation and the a effect is 
induced by azimuthally drifting Rossby waves. 
Figures [2] and [3] show the dynamo solution at E = 10~ 3 ' 5 , Ro = —1 and Rm c = 2650. Figure 



2(a) shows the axisymmetric flow and field in the meridional plane and figure 2(b) the flow 



and field distributions in the equatorial plane. The angular velocity and meridional circulation 



in figure 2(a) exhibit the Stewartson layer and the process of Ekman pumping. The flow in the 
equatorial plane indicates that the Stewartson layer has already been destabilised at Ro = — 1 
and the m = 1 mode of an azimuthally drifting wave has emerged, which is consistent with the 
linear stability calculations for spherical Couette flow (right panel in figured]). The poloidal 
magnetic field exhibits a dipolar symmetry about the equator and in the equatorial plane the 
field exhibits an m = 1 mode. 

To better understand this dynamo we diagnose some percentages of kinetic and magnetic 
energies. Figure [3] shows the spectrum of kinetic and magnetic energies. The percentage of 
toroidal kinetic energy is 98.6% (Tabled]) and the percentages of kinetic energies in m = and 
m = 1 modes are, respectively, 98.1% and 1.8% (upper right panel in figure [3]). Clearly the 
axisymmetric azimuthal flow is dominant and the differential rotation of this flow shears the 
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(a) Contours of axisymmetric flow and field in the meridional plane. From left to right, 
angular velocity U^/ (r said), meridional circulation <!/ satisfying V X (^(^/(r sin#)) = 
(U r ,Ug), toroidal field and lines of poloidal field \ satisfying V X (\4>/{r sind)) = 
(B r ,Bg). In the panels of angular velocity and toroidal field, solid lines represent west- 
east direction and dashed lines east- west. In the panels of meridional circulation and 
poloidal field, solid lines represent anti-clockwise direction and dashed lines clockwise. 
Contour lines are uniformly spaced with 10 levels. 




(b) Contours of flow and field in the equatorial plane. From left to right, radial (U r and 
B r ), colatitude (Ug and Bg) and azimuthal (Uj, and Bj,) components of flow (top row) 
and field (bottom row). Solid lines represent positive value and dashed lines negative. 
Contour lines are uniformly spaced with 10 levels. 



Figure 2. E = 10~ 35 , Ro = -1 and Rm c = 2650. 

poloidal field to create a strong toroidal field, i.e. the uj effect. This assertion is supported by 
the dominance of toroidal magnetic energy, namely that the percentage of toroidal magnetic 
energy is 97.8% (Table [T]). The percentages of magnetic energies in the m = and m = 1 
modes are, respectively, 81.4% and 15.6% (lower right panel in figure[3|). So this dynamo is an 
a-bj dynamo, and the a effect is induced by Rossby waves which mainly consist of the m = 1 
mode (here we address again that the a effect does not connect to the mean field theory). A 
weak flow in the m = 1 mode (1.8%) can generate a relatively stronger poloidal field (2.2%), 
so this Rossby- wave- induced a effect is efficient. The efficiency may be interpreted through 
the location of the a and uj effects. Figure [2(b)| shows that Rossby waves are concentrated in 
the vicinity of the inner sphere where the Stewartson layer is located and where there is a 
strong differential rotation to induce the uj effect. Therefore the a effect arising from Rossby 
waves and the uj effect arising from differential rotation work at the same location, which 
strengthens dynamo action. The efficiency resulting from location is discussed in detail by 
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Figure 3. Energy spectrum of flow and field for the solution in figure [2] The top row shows the kinetic energy and the 
bottom row the magnetic energy. The left column shows the energy as a function of the degree I of spherical harmonics 
and the right column the energy as a function of the order m. The large figures are plotted using a linear scale and 
the inset figures in the left column embedded in the large figures are plotted using a logarithmic scale. The solid lines 
represent the total energy, the dashed lines the toroidal energy and the dot-dashed lines the poloidal energy. 



Gubbins and Gibbons (120091 1. 

We then change the sign of Ro to be positive, i.e. prograde flow. Figures |4] and [5] show the 
dynamo solution at E = 10 -3 ' 5 , Ro = +1 and Rm c = 5000. The critical Rm for this prograde 
flow is greater than that for the retrograde flow. This asymmetry of the onset of dynamo 
action arises from the asymmetry between the prograde and retrograde flows, as we shall now 
illustrate. 

Firstly, in the axisymmetric part of the two flows, we can observe that the Ekman layer of 
the Ro = +1 flow is thinner than that of the Ro = — 1 flow (the first panels in figures 2 (a 



and 4(a)). When \Ro\ is very small, the thickness of the Ekman layer 5 depends only on E, 
i.e. S ~ \fE. This classic power law is derived from the balance among the viscous force, the 
Coriolis force and the pressure gradient with the neglect of the inertial force in the laminar 
Ekman layer (flows in our calculations are supercritical but still laminar). But when \Ro\ is 
finite the inertial force is not negligible and so 5 also depends on Ro. The reason that 5 of 
the Ro = +1 flow is thinner than that of the Ro = — 1 flow is because of the average absolute 
rotation rate. The average absolute rotation rate of the Ro = +1 flow is between +1 and 
+2 whereas that of the Ro = —1 flow is between +1 and 0. A stronger rotation rate induces 
a thinner Ekman layer. Moreover, the Ekman layer of the Ro = +1 flow is thickened near 
the equator whereas that of the Ro = — 1 flow is squashed. This is because of the opposite 



directions of the meridional circulation (the second panels in figures 2(a) and 4(a) ). Near the 
equator there is a radially outward jet in the Ro = +1 flow but a radially inward jet in the 
Ro = — 1 flow. 

Secondly, according to figured! retrograde flow is more stable than prograde flow such that 
the Ro = +1 flow is more supercritical than the Ro = — 1 flow (Ro c = +0.5 for prograde 
flow and Ro c = —0.8 for retrograde flow). This can be seen from the two flow patterns. In the 
meridional plane both the angular velocity and the meridional circulation of the Ro = +1 flow 
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(a) As in figure |2(a)| 





(b) As in figure [2(b)] 
Figure 4. E = 10~ 3 - 5 , Ro = +1 and Rm c = 5000. 



have more complex structures than those of the Ro = — 1 flow. In the equatorial plane the 
Ro = +1 flow exhibits an m = 3 mode which is consistent with the linear stability calculations 
for E = 10~ 3 ' 5 (left panel in figure [1]). This mode occupies 14.7% of the kinetic energy (upper 
right panel in figure [5]) which is much larger than the 1.8% of the m = 1 mode in the Ro = — 1 
flow. Moreover, the spiral structure of the Rossby waves for prograde flow can be clearly 
perceived to spread radially outward (top row in figure 4(b) [ ) whereas the Rossby waves for 
retrograde flow are in the vicinity of the inner sphere (top row in figure 2(b) ). 

All these differences between prograde and retrograde flows in respect of both axisymmetric 
and nonaxisymmetric parts can be mathematically interpreted through the Coriolis force. The 
Coriolis force is a first order term in the N-S equation such that it will bring opposite effects 
by changing the sign of the flow. Therefore it is not surprising that prograde flow requires 
a larger Rm c than retrograde flow, because in retrograde flow the a effect induced by the 
Rossby waves works near the tangent cylinder, where the u effect induced by differential 
rotation works, but in prograde flow the Rossby waves spread radially outward and therefore 
the a and u effects work at different locations, which weakens dynamo action. 

The magnetic field i n thi s solution has a rich spectrum governed by the selection rules of 
Bull ard and Gellman (I1954T ). Since the flow is m = 3 no harmonics of the form m = 3i (i = 
1,2,.. .) can be generated, which is shown in the lower right panel in figure O This is more 
complex than the gently decaying spectrum for Ro = — 1 shown in the lower right panel of 
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Figure 5. As in figure [3] but for the solution in figure [4] The I = 2 contribution of flow is identically zero and thus is 
excluded from the log-log plot (inset figure in the upper left panel). 

figure O which is dominated by m = and m = 1. The field has a quadrupolar symmetry 
about the equator. The equatorial symmetry depends on the parity of (I — m) where I and 
m are the degree and order of the spherical harmonic components of the poloidal field. For 
the axisymmetric field, m = and then the equatorial symmetry is determined by Z. In the 
Ro = — 1 dynamo solution the dominant I is odd, so the axisymmetric field is predominantly 
dipolar. However, in the Ro = +1 dynamo solution the dominant I is even, so the axisymmetric 
field is predominantly quadrupolar. 

Next we reduce the Ekman number to 10 -4 and we firstly test retrograde flow. We try 
to use the flow at Ro = — 1 to produce a dynamo, but unfortunately we fail to do so. Until 
Rm = 10 we cannot find any dynamo. This is con sistent with the nonlinear spherical Couette 
dynamo calculations ( Guervillv and Car din 2010l ) in which dynamo action in retrograde flow 
at E = 10 -4 does not occur until \Ro\ is increased b eyond 0.5 which is tr a nslate d to be 
1.4 in our definition. Here we should notice that E in iGuervillv and Cardinl (|20ld ) and in 
our definition differs by a factor (1 — rj/r D ) 2 (equation fl5H)), and so this comparison is not 
strictly correct, but it gives us a qualitative result that retrograde flow at low E does not 
favour dynamo action. We may increase \Ro\ to be greater in a more supercritical regime, say, 
Ro < — 1 in which the absolute rotation rate of inner sphere is in the opposite direction, but 
our computational power is not capable to do that. 

We now move to prograde flow. We choose Ro to be +0.3 which is slightly above the neutral 
stability curve (left panel in figured]) and we find dynamos. Figures [6] and [7] show the dynamo 
solution at E = 10~ 4 , Ro = +0.3 and Rm c = 6500. Both the flow and field exhibit a more 
columnar structure than E = 10~ 3 ' 5 (figures 6(a) and 4(a) ) due to the stronger Coriolis effect 
(the Taylor-Proudman theorem) . The flow exhibits a nice spiral structure and an m = 3 Rossby 
wave spreads radially outward. The m = 2 and m = 1 modes of the magnetic field occupy, 
respectively, 40.0% and 33.9% of the total magnetic energy (lower right panel in figure [7]) and 
the field has a quadrupolar symmetry because the dominant I of the axisymmetric poloidal 
field is even. 
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(a) As in figure |2(a)| 




(b) As in figure |2(b) | 



Figure 6. E = 10" 4 , Ro = +0.3 and Rm c = 6500. 

The reason that at E = 10 -4 a retrograde flow (Ro = — 1) cannot easily produce a dynamo 
but a prograde flow (Ro = +0.3) can is because of the flow supercriticality. The asymmetry of 
linear stabilities between progarde and retrograde flows increases as the global rotation rate 
increases (E decreases, see figure [I]). According to figure [H at E = 10 -4 the magnitude of the 
critical Ro of retrograde flow is 0.4 (Ro c = —0.4) but that of prograde flow is 0.2 (Ro c = +0.2). 
The former being twice as large as the latter, we may postulate that the flow supercriticality 
at Ro = +0.3 is already strong enough to generate a dynamo but the flow supercriticality 
at Ro = — 1 is not yet sufficient. For example, the percentage of nonaxisymmetric modes in 
kinetic energy is 6.68% for the Ro = +0.3 flow whereas it is 6.56% for the Ro = — 1 flow. 

To verify this point that the strong flow supercriticality favours dynamo action, we continue 
to increase Ro from +0.3 to +1. Figures [8] and [9] show the dynamo solution at E = 10 -4 , 
Ro = +1 and Rm c = 2000. That Rm c at Ro = +1 is much lower than Rm c at Ro = +0.3 
supports the point that the stronger flow supercriticality is better for dynamo action. On the 
other hand, comparison between Rm c = 2000 for (E = 10~ 4 ,i?o = +1) and Rm c = 5000 for 
(E = 10~ 3 ' 5 , Ro = +1) also supports this point. Ro c for E = 10~ 4 is +0.2 whereas Ro c for 
E = 10~ 3 ' 5 is +0.5 (figured]), and then Ro = +1 results in stronger flow supercriticality for 
the E = 10 -4 flow than for the E = 10 -3 ' 5 flow, such that Rm c of the E = 10 -4 flow is lower 
than Rm c of the E = 10~ 3 - 5 flow. 

Although the Ro = +1 flow also exhibits an m = 3 mode (top row in figure [8(b)] and upper 
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Figure 7. As in figure [3] but for the solution in figure [6] The I = 2 contribution of flow is identically zero and thus is 
excluded from the log-log plot (inset figure in the upper left panel). 

right panel in figure [9]) as does the Ro = +0.3 flow, it has a more complex structure (top rows 
in figures [8(b)1 and [6(b)|) arising from the stronger nonlinear inertial force. Accordingly the 
magnetic field has a more complex structure (bottom rows in figures [8(b)1 and [6(b)] ). Therefore 
we can conclude that the stronger flow supercriticality which leads to more complex struc- 
tures in the flow and field facilitates the onset of dynamo action because of the more efficient 
interaction between different modes, namely the induction term V x (17 x B) in the magnetic 
induction equation ([2]). Although some simultations in turbulent regime showed that an in- 
crease of turbulence increase d the dynamo thres hold due to the enhancement of the magnetic 
diffusivity by the turbulence ( Bavliss et al. 20071 ). our calculations in this supercritical regime 
are far away from the turbulent regime. 

Following our calculations for different E and Ro (Ro in respect of both sign and magni- 
tude), we implement a perfectly electrically conducting boundary condition to investigate the 
influence of magnetic boundary conditions on dynamo action. Figure [10] shows the dynamo 
solution at E = 10~ 3 ' 5 , Ro = — 1 and Rm c = 2450 with both spheres being perfectly con- 
ducting. We show only the axisymmetric flow and field distributions in the meridional plane. 



The flow/field distribution in the equatorial plane is identical/similar to figure 2(b) and so 
is not shown again. Compared with Rm c = 2650 for insulating boundaries in figure [2] the 
conducting boundaries favour dynamo action, at least in this parameter regime of E = 10~ 3 ' 5 
and Ro = — 1. This may be interpreted with the argument that in the case of conducting 
boundaries the poloidal field lines are trapped in the spherical shell and cannot penetrate out 
(the fourth panel in figure [TO]) such that no poloidal magnetic energy escapes, whereas the 
poloidal field lines with insulating boundaries can penetrate out (the fourth panel in figure 
2(a) ) such that some poloidal magnetic energy is lost. 

In addition to the magnetic boundary condition, we perform calculations using a different 
aspect ratio rj/r D = 1/2 to investigate the influence of the aspect ratio on dynamo action. 
Figure [TT] shows the dynamo solution at E = 10~ 3 ' 5 , Ro = —1 and Rm c = 1250 with 
r i/ r o = 1/2- The flow exhibits an m = 1 mode and the field has three comparable modes 
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(a) As in figure |2(a)| 




(b) As in figure |2(b) | 
Figure 8. E = 10~ 4 , Ro = +1 and Rm c = 2000. 



m = 0, 1 and 2. The critical magnetic Reynolds number Rm c = 1250 for r,/r D = 1/2 is lower 
than Rm c = 2650 for rj/r D = 1/3. This dynamo action is more efficient. One reason for this 
is that the total kinetic energy for rj/r D = 1/2 is stronger than that for rj/r Q = 1/3. In 
our dimensionless equations the spherical gap L is unity, and so the radii of the spheres are 
r a = 1.5 and rj = 0.5 for rj/r G = 1/3, and r Q = 2 and r« = 1 for r,/r D = 1/2. According 
to equation ([6b), the inner boundary velocity for rj/r D = 1/2 is twice that for r,i/r = 1/3 
such that a stronger flow is driven with ri/r = 1/2. Note that our current definition of the 
magnetic Reynolds number Rm is based on the spherical shell width (equation j3b)). When 
one translates this into a new magnetic Reynolds number Rm' based on the outer sphere 
radius given by 



Rm' 



\A(2\ 



:ila) 



the following interrelation holds: 



Rm = Rm 



1 - — 



(lib) 



Based on the outer shell magnetic Reynolds number, the critical Rm' are 5960 and 5000 
respectively for the cases ri/r = 1/3 and rj/r G = 1/2. Thus one can argue that the small 
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Figure 9. As in figure [3] but for the solution in figure [8] The I = 2 contribution of flow is identically zero and thus is 
excluded from the log-log plot (inset figure in the upper left panel). 







Figure 10. E = 10 , Ro = —1 and Rm c = 2450, and perfectly electrically conducting boundaries. 



change in aspect ratio has led to a more efficient dynamo and a lower dynamo threshold. The 
reason for this increased efficiency is the strength of the a effect: in this dynamo solution 
with ri/r = 1/2, the percentages of kinetic energies in the m = and m = 1 modes are, 
respectively, 96.1% and 3.5%. Compared with the percentage of 1.8% for the m = 1 mode in 
the dynamo solution with ri/r = 1/3, the percentage of Rossby waves for the ri/r a = 1/2 
case is doubled such that the strength of the a effect is doubled . This point can be also 
supported by Re w and Re a in Tabled! In Hollerbach et al\ ( 20061 ) various aspect ratios for 
spherical Couette flow have been studied and the aspect ratio indeed alters the pattern of 
Rossby waves. This can be interpreted as a result of the /3 effect, in which the curvature of 
the boundary induces Rossby waves. f3 is defined as f3 = l/h(dh/ds) where h is the height of 
the fluid column at cylindrical radius s. Clearly the value of f3 is related to the aspect ratio. 
In the shallow wa ter theory of a thin spherical gap, (5 appears in the dispersion relationship 
for Rossby waves (lPedloskylll987l ). and although in the geometry of a wide spherical gap there 
is no explicit dispersion relationship, (3 is still essential to Rossby waves. Therefore it is not 
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surprising that rj/r Q = 1/2 is better for dynamo action because both the total flow and the 
a effect are stronger. 




(a) As in figure |2(a)| 









(b) As in figure [2(b)] 
Figure 11. E = 10 -3 ' 5 , Ro = —1 and Rm c = 1250, and aspect ratio ri/r = 1/2. 



For a better understanding of the "dynamo window" in the nonlinear spherical Couette 
dynamo ( Guervillv and Cardin 2010), we test the corresponding kinematic spherical Couette 
dynamo. At first we translate the paramters of "dynamo window" to our definition and they 
are E = 10" 2 - 6 , Ro = -4.3 and Rm c ~ 360. Then we calculate with E = lO" 2 6 , Ro = -4.3 
and Rm = 400. However, we cannot find a dynamo. Therefore, we may speculate that this 
low Rm "dynamo window" is mathematically a subcritical dynamo in the nonlinear regime. 

To end this section we investigate the direction of magnetic dipole axis in all the six dynamo 
solutions. We compare between the contributions of two modes (I = 1, m = 0) and (I = 1, m = 
1) of poloidal component, i.e. h\o and h\\ (equations ([8]) and ([9])), to the radial field at the outer 
boundary. The former represents the axial dipole whereas the latter the equatorial dipole. It is 
interesting that the magnetic dipole is almost axial in the retrograde flows whereas it is almost 
equatorial in the prograde flows. It is not surprising because the magnetic field generated by 
the retrograde flows has very small portion of energy in the m = 1 mode (lower right panel in 
figure [3]) whereas the field generated by the prograde flows has very small portion of energy 
in the m = mode (lower right panels in figures 02 [7] and [9]) . But we do not understand the 
physics why the direction of Ail is crucial to the direction of magnetic dipole axis. 
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4. Conclusions 

In this work we have calculated examples of kinematic dynamos driven by spherical Couette 
flow. We draw the following conclusions. 

2D axisymmetric flows at the Ekman number regime in our calculations cannot generate 
dynamo action whereas 3D flows can. This indicates that azimuthally drifting Rossby waves 
are crucial to dynamo action. The spherical Couette dynamos we have studied are a-u dy- 
namos in which the oj effect is induced by differential rotation in the Stewartson layer and the 
a effect is induced by helical Rossby waves arising from the destabilisation of the Stewartson 
layer. 

At a slow global rotation rate (E = 10~ 3 ' 5 ), retrograde flow favours dynamo action but 
prograde flow does not, because in the former the a and oj effects work at the same location 
near the tangent cylinder but in the latter they two work at different locations. At a rapid 
global rotation rate (E = 10 -4 ), prograde flow favours dynamo action but retrograde flow 
does not, because in the former the flow supercriticality arising from the nonlinear inertial 
force is strong enough to generate a dynamo whereas in the latter it is not sufficient. With 
a higher global rotation rate this asymmetry in the flow supercriticality is more prominent. 
Comparison between Ro = +0.3 and +1 at the same E = 10 -4 and comparison between 
E = 10 -3 ' 5 and 10~ 4 at the same Ro = +1 both suggest that stronger flow supercriticality 
(leading to more complex structures of flow and field) facilitates the onset of dynamo action. 

A conducting magnetic boundary condition and a larger aspect ratio both favour dynamo 
action. Through our calculation of kinematic spherical Couette dynamo we may speculate 
that the low critical Rm "dynamo window" in the nonlinear spherical Couette dynamo could 
be mathematically interpreted as a subcritical dynamo. The magnetic dipole is almost axial in 
the retrograde flows and almost equatorial in the prograde flows. We also hope that this work 
may be of help to the o ngoing experiments of spherical Couette dynamo, e.g. the dynamo 
experiment in Maryland (IPefflev et aZJl200flh . 
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